First-Principles Study of Electron Linewidths in Graphene 
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We present first-principles calculations of the linewidths of low-energy quasiparticles in n-doped 
graphene arising from both the electron-electron and the electron-phonon interactions. The contri- 
bution to the electron linewidth arising from the electron-electron interactions vary significantly with 
wavevector at fixed energy; in contrast, the electron-phonon contribution is virtually wavevector- 
independent. These two contributions are comparable in magnitude at a binding energy of ~ 0.2 eV, 
corresponding to the optical phonon energy. The calculated linewidths, with both electron-electron 
and electron-phonon interactions included, explain to a large extent the linewidths seen in recent 
photoemission experiments. 



Graphene [H, Si, a single layer of carbon atoms in 
a hexagonal honeycomb structure, is a unique system 
whose carrier dynamics can be described by a massless 
Dirac equation [4!]. Within the quasiparticle picture, car- 
riers in graphene exhibit a linear energy dispersion rela- 
tion and chiral behavior resulting in a half-integer quan- 
tum Hall effect Q absence of backscattering @], 
Klein tunneling [7|, and novel phenomena such as elec- 
tron supercoUimation in superlattices [§, [l^ . 

Graphene is considered a promising candidate for elec- 
tronic and spintronic devices [llj . For these applications 
it is important to understand the effects of many-body in- 
teractions on carrier dynamics. In particular, the scatter- 
ing rate of charge carriers, manifested in their linewidths, 
affects the transport properties of actual devices. 

The scattering of charge carriers in solids can 
arise from several different mechanisms, among which 
electron-hole pair generation, electron-plasmon interac- 
tion, and electron-phonon (e-ph) interaction are gener- 
ally important. Scattering by impurities, defects and in- 
teractions with the substrate also affects the carrier dy- 
namics. The contribution to the electron linewidths aris- 
ing from the e-ph interaction has been studied with first- 
principles calculations , Ts'l and through the use of an- 
alytical and numerical calculations based on the massless 
Dirac equation [iJi EBl • The linewidth contribution origi- 
nating from electron-electron (e-e) interactions, which in- 
cludes both the electron- hole pair generation process and 
the electron-plasmon interaction, has only been studied 
within the massless Dirac equation formalism [TgI. [TtI. [isj . 

A recent angle-resolved photoemission experiment on 
n-doped graphene epitaxially grown on silicon carbide 
(SiC) [19|] has stimulated experimental [H [U, and 
theoretical [H, Q [H, [H, ElUl studies on this topic. In 
Ref. Tl, the width of the momentum distribution curve 
(MDC) from photoemission data is presented. The MDC 
of the graphene photoemission spectra is observed to 
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resemble a simple Lorentzian whose width may be in- 
terpreted to be directly proportional to the scattering 
rate [H. 

We draw the attention to the well-known controversy 
in the different interpretations of the angle-resolved pho- 
toemission spectra of graphene. It is claimed in Ref. [l^ 
that the spectral features can entirely be understood from 
many-body effects, including both e-e and e-ph inter- 
actions, in graphene. On the other hand, in Ref. [20l . 
it is argued that many of those features are dominated 
by an energy gap of 0.2^0.3 eV, which opens up at 
the Dirac point energy (Sd) because of interactions be- 
tween graphene and the reconstructed surface of SiC. 
This important problem in understanding the quasiparti- 
cle spectra of graphene (which also have implications in 
graphene-based electronics applications) has led to nu- 
merous additional experiments directly or indirectly ad- 
dressing this discrepancy [H, [13, [H, [11, [13, HI . On the 
theoretical side, several density functional theory calcu- 
lations on the effect of substrates without considering 
many-body effects, along the line of Ref. '2^, have been 
performed 29, 30, 31]. On the other hand, first-principles 
calculations on the effects of both e-e and e-ph interac- 
tions, along the line of Ref. [l^, have been lacking up to 
now. 

In this paper, to fill in this missing part, we present ab 
initio calculations of the electron linewidth in n-doped 
graphene arising from e-e interactions employing the GW 
approximation |32l . [33l . [3^ . In addition, we calculate the 
electron linewidth originating from the e-ph interaction 
following the method in Rcfs. [T^ [35l and [36l. Combining 
both contributions, we provide a comprehensive view of 
the scattering rate originating from many-body effects. 
Our calculation indicates that the linewidth arising from 
e-e interactions is highly anisotropic. This is in contrast 
to the insensitivity to wavevector of the phonon-induced 
electron linewidth shown in Ref. [H. The calculated 
linewidth arising from e-e interaction becomes compa- 
rable to that arising from e-ph interaction at a binding 
energy of ^ 0.2 eV (i.e., the optical phonon energy). 
The combination of the two contributions accounts for 
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most of the measured linewidth over the eV ~ 2.5 eV 
binding energy range. 

The electronic eigenstates |nk) of graphene are ob- 
tained with ab initio pseudopotential density-functional 
calculation s | 37I| in the local density approximation 
(LDA) [sl, l39f in a supercell geometry. Electronic wave- 
functions in a 72 X 72 X 1 k-point grid are expanded in 
a plane-waves basis (40| with a kinetic energy cutoff of 
60 Ry. The core-valence interaction is treated by means 
of norm-conserving pseudopotentials f4l| . Graphene lay- 
ers between adjacent supercells are separated by 8.0 A 
and the Coulomb interaction is truncated to prevent spu- 
rious interaction between periodic replicas [42| . Increas- 
ing the interlayer distance to 16.0 A makes virtually no 
difference in the calculated electron self-energy. Doped 
graphene is modeled by an extra electron density with a 
neutralizing background. 

We calculate the imaginary part of the electron self- 
energy induced by the e-e interaction within the GFK ap- 
proximation [33, 113] . The frequency dependent dielectric 
matrices eG,G'(q, t^) are calculated within the random 
phase approximation using the LDA wavefunctions on a 
regular grid of iv with spacing Aiv =0.125 eV [i^, and 
a linear interpolation is performed to obtain the dielec- 
tric matrices for energies in between the grid points. In 
the calculation of the polarizability, for numerical con- 
vergence, an imaginary component of magnitude Auj of 
0.125 eV as above is introduced in the energy denom- 
inator. Convergence tests showed that the dimension 
of ec.G' may be truncated at a kinetic energy cutoff 
of ?i^G^/2m =12 Ry. To take into account the screen- 
ing of the SiC substrate, we have renormalized the bare 
Coulomb interaction by an effective background dielec- 
tric constant of £b = (1 + esic)/2 = 3.8 [3 113, Elil , 
where esio(= 6.6) is the optical dielectric constant of 
SiC 0,1^. 

Figure [1] shows the calculated imaginary part 
IniE,''^j^''(£„k) = (j^kl Inil]''~''(r,r',£„k) |"-k) of the elec- 
tron self-energy arising from the e-e interaction with uj set 
at the LDA eigenvalue e„k- The Fermi level Ep (= 0) is 
taken to be 1 eV above Eu- In Fig.[lja), Im S,'^jk'^(enk) for 
graphene without including substrate screening, appro- 
priate for suspended graphene [H, \^ , is plotted along 
the KF direction. Generally, the self-energy increases 
with increasing |£„k| as measured from Ep. A notable 
feature is the peak around £„k = —1.5 eV. To find the 
origin of this peak, we have decomposed the total electron 
self-energy into the contributions arising from transitions 
into the upper linear bands (above Eu) and the lower 
linear bands (below -Ed)- The former involves electron- 
plasmon interaction . The peak structure comes from 
scattering processes of electrons into the upper linear 
bands, whereas those scattering processes into the lower 
linear bands result in a monotonic increase in the electron 
linewidth. When the background dielectric constant £b 
is changed from 1 to 3.8 [Fig. [Ub)], the position of this 
peak shifts toward lower-binding energy by ~ 0.3 eV, 
reflecting a decrease of the plasmon energy in graphene 



> 

°i. 0.3 



r -\ 
-\ ^ *7 K ^ 

'v y \ 




- (a) 

^£b=1.0^ 


>-<,^ 


Ed^ 

( r 


- (b) 

^£b=3.8 : 

"\M / 


< r >K \ 




- (c) 




( r 


- (d) 

^8b=3.8 

'\M 
>K 



-3 -2-10 1 -2-10 1 
Electron energy 8„i^ (eV) Electron energy e„|^ (eV) 




0.5 1 
h vq q (eV) 



FIG. 1: (color online), (a)-(d): Calculated imaginary part 
of the electron self-energy arising from the e-e interaction, 
Im E^^'^(£„k), versus the LDA energy e„k (solid black lines) in 
n-doped graphene. The Dirac point energy Eu is 1.0 eV below 
the Fermi level. The contributions to ImE^^'^(e„k) coming 
from electronic transitions to the upper linear bands and to 
the lower linear bands are shown as dashed red lines and dash- 
dotted blue lines, respectively. The self-energy is evaluated 
along the reciprocal space segments shown in the insets, (a) 
and (c) are results for suspended graphene with a background 
dielectric constant of eb = 1-0, whereas (b) and (d) are results 
for graphene with a background dielectric constant of = 
(1 + esic)/2 = 3.8. The Fermi level and Eu are indicated 
by vertical lines, (e): Calculated plasmon energy dispersion 

relation given by eG=o,G'=o[q, <^oo(ci)] = 0' versus 

?!.uo|q| along the FM direction. The solid lines are guides to 
the eye and the dashed line corresponds to i-o[q) — hvoq. 

[Fig. (He)] m, m. The height of the peak is further 
suppressed. At low energy (|£nk| ^ 1-0 ^V); the imag- 
inary part of the self-energy is however not sensitive to 
the choice of £b. 

Comparing Figs.[Ha) and[2b) with Figs.[T]Jc) and[Ud) 
shows that the electron self-energy arising from the e- 
e interaction calculated along the KM direction is very 
different from that along the KF direction. Below 
— 1.5 eV, Iml]^^'^(£„k) along K— >M decreases with in- 
creasing |£nk|, and it almost vanishes at the M point. 
This strong k anisotropy in the e-e contribution to the 
imaginary part of the self-energy is a band structure ef- 
fect, and is absent in calculations based on the massless 
Dirac equation. This behavior is in contrast with the 
wavevector insensitivity of the phonon-induced electron 
self-energy flS\ (Fig. The calculated real part [l^l 
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FIG. 2: (color online). Calculated Im E„k(£nk) versus the 
LDA energy eigenvalue e„k in n-doped graphene (i?D = 
— 1.0 eV) on a model substrate (sb = 3.8). The total self- 
energy, the self-energy arising from the e-e interaction, and 
that arising from the e-ph interaction are shown in solid black, 
dashed red and dash-dotted blue lines, respectively. The self- 
energy is evaluated along the reciprocal space segments shown 
in the insets. 



and the imaginary part (SJ] of the electron self-energy 
in bulk graphite arising from the e-e interaction are also 
anisotropic, in line with the present findings. 

Figures HJa) and [DJb) show the electron self-energy in 
n-doped graphene [E-q = — 1.0 eV) on a substrate (model 
with Eb = 3.8) arising both from the e-e and the e-ph 
interaction. The Im S„k(enk) along the two different di- 
rections KF and KM are qualitatively different at high 
binding energy. This anisotropy is due to the e-e in- 
teraction, and not the e-ph interaction [l3| . It is noted 
that the total linewidth along the KM direction is almost 
constant for binding energies in the range 1.7 to 3.5 eV. 
These anisotropic features should be observable in pho- 
toemission experiments. 

The e-e and the e-ph interactions give comparable con- 
tributions to the imaginary part of the electron self- 
energy, especially within a few tenths of an eV from 
the Fermi level (Fig. [2]). This behavior is peculiar 
to graphene. In most metals the e-ph contribution 
to the electron self-energy near E-p is generally domi- 
nant over the e-e contribution at e nerg ies comparable to 
the relevant phonon energy scale [5ll |. Similarly large 
Im S„k(£nk) due to e-e interactions are obtained in the 
Dirac Hamiltonian calculations in Refs. [l^ and [13 if the 
same background dielectric constant is used. Because 
of this peculiar aspects of graphene, an e-ph coupling 
strength A extracted from measured data could be over- 
estimated if the e-e interaction is neglected. This may 
explain why the e-ph coupling strength A extracted from 
photoemission spectra [2l| is larger than the theoretical 
calculations [T^ together with the effects of bare 
band curvature [l3| and dopants. 

We now compare the imaginary part of the elec- 
tron self-energy obtained from our calculation with the 
MDC width obtained from measured photoemission spec- 
tra p^. For a linear bare band energy dispersion, the 
spectral function at a fixed energy u; is a Lorentzian 
as a function of the wavevector measured from the K 
point [l^. Thus, the width of the MDC Afc at energy lo = 
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FIG. 3: (color online). MDC width versus binding energy 
in n-doped graphene {Eo = —1.0 eV). Calculated quanti- 
ties for suspended graphene (cb = 1-0) and for graphene on a 
model substrate (sh ~ 3.8) are shown in dash-dotted blue and 
dashed red lines, respectively. The experimental result mea- 
sured for sample corresponding to the highest level of doping 
in Fig. 3 of Ref. JJ2 are shown as the solid black line [s^l ■ Both 
the experimental and the calculated results are along the KM 
and the KF direction of the Brillouin zone when the electron 
energy is above and below -Ed, respectively. 



e„k can be identified as Afc(£„k) = 2ImI]„k(£nk)/?it'o 
where vq is the LDA band velocity of low-energy charge 
carriers in graphene [HI, [l^ . (For the n-doped graphene 
with Eu ~ —1.0 cV, the bare band dispersion is, to a 
good approximation, linear in the energy range consid- 
ered in Fig. O) 

Figure [3] shows the calculated MDC width for sus- 
pended graphene (sb = 1-0) and for our model of 
graphene on SiC (eb = 3.8). The substrate screen- 
ing affects the position and the strength of the peak 
arising from the electron-plasmon interaction, while the 
low-energy part is insensitive to the dielectric screen- 
ing from the substrate. The calculated MDC width 
for graphene when substrate screening is accounted for 
is in agreement with the experimental data of Ref. [l^ 
throughout the whole energy window shown in Fig. [31 
However, the experimentally measured MDC width in a 
0.4 eV energy window around i?D (=-1.0 cV) is larger 
than that from our calculation. This enhanced linewidth 
may possibly arise from the gap which opens up at i?D 
and midgap states originating from the interactions be- 
tween graphene and SiC substrate with a carbon buffer 
layer [20, |2l H HI . 

In conclusion, we have studied the electron linewidths 
of n-doped graphene including both the e-e and the e- 
ph interaction contributions, using first-principles calcu- 
lations. The imaginary part of the electron self-energy 
arising from the e-e interaction is strongly anisotropic 
in k-space. We have shown that for graphene, unlike in 
conventional metals, the e-e contribution is comparable 
to the e-ph contribution at low binding-energy. Our cal- 
culation explains most of the scattering rate observed in 
a recent photoemission experiment [19|; however, near 
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the Dirac point energy, the calculated scattering rate is 
smaller than the measured one, suggesting the possibil- 
ity of band gap opening and midgap states. These re- 
sults contribute to the resolution of the important con- 
troversy introduced earlier in this paper and encourages 
further theoretical studies including both many-body in- 
teractions and substrate effects at an atomistic level. 
More generally, our first-principles calculations convinc- 
ingly demonstrate that multiple many-body interactions 
ought to be considered on the same footing in order to 
achieve a quantitative and comprehensive interpretation 
of high-resolution angle-resolved photoemission spectra. 
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